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Pl,- ABSTRACT 

O ' Galaxy harassment is an important mechanism for the morphological evolu- 

& ' 

' tion of galaxies in clusters. The spiral galaxy NGC 4254 in the Virgo cluster 

is believed to be a harassed galaxy. We have analyzed the power spectrum of 
HI emission fluctuations from NGC 4254 to investigate whether it carries any 

> 

qq ■ imprint of galaxy harassment. The power spectrum, as determined using the 



16 central channels which contain most of the HI emission, is found to be well 
fitted by a power law P(U) = AU a with a = — 1.7 ± 0.2 at length-scales 
. 1.7 kpc to 8.4 kpc. This is similar to other normal spiral galaxies which have a 

slope of ~ —1.5 and is interpreted as arising from two dimensional turbulence 
at length-scales larger than the galaxy's scale-height. NGC 4254 is hence yet 
another example of a spiral galaxy that exhibits scale-invariant density fluctu- 
ations out to length-scales comparable to the diameter of the HI disk. While 
a large variety of possible energy sources like proto-stellar winds, supernovae, 
shocks, etc. have been proposed to produce turbulence, it is still to be seen 
whether these are effective on length-scales comparable to that of the entire 
HI disk. On separately analyzing the HI power spectrum in different parts of 
NGC 4254, we find that the outer parts have a different slope (a = — 2.0±0.3) 
compared to the central part of the galaxy (a = —1.5 ± 0.2). Such a change 
in slope is not seen in other, undisturbed galaxies. We suggest that, in addi- 
tion to changing the overall morphology, galaxy harassment also effects the 
fine scale structure of the ISM, causing the power spectrum to have a steeper 
slope in the outer parts. 



Key words: physical data and process: turbulence-galaxy:disk-galaxies:ISM 



2 Prasun Dutta, Ayesha Begum, Somnath Bharadwaj and Jayaram N. Chengalur 
1 INTRODUCTION 



Galaxy harassment (frequent high speed galaxy encounters, iMoore et al.l Il996l ) is believed 
to be an important process in driving the morphological transformation of spiral galaxies 
to ellipticals inside clusters. Typically, the first encounters convert a normal spiral galaxy 
to a disturbed spiral with dramatic features drawn out from the dynamically cold gas. 
he spiral galaxy N GC 4254, located in the nearby Virgo cluster, is found to have a tail 



( IMinchin et al.ll2005l ) with neutral hydrogen (HI) mass 2.2 x 10 8 M Q within a distance of 



120 kpc from the galaxy. This gaseous tail, without an y stellar counterpar t, is believed 



to have been produced by an act of galaxy harassment (lHaynes et al.l 120071 ). Each act of 
harassment has the potential to induce a burst of star formation and to change the internal 
properties of the galaxy, including the properties of the inter stellar medium (ISM). Here we 
study the effect of the harrasment on the fine scale structure of the ISM. We use the power 
spectrum of HI intensity fluctuations to quantify the fine scale structure. 



Power spect rum analysis of the HI 2 



cm intensity fluctuation has been widely used to 



the ISM ( 


Crovisier & Dickev 


1983; 


Green 


1993; 


Stutzki et al. 


1998: 


Stanimirovic et al. 




Deshpande et al. 


2000; 


Elmeexeen et al. 


2001 


Elmeereen & Scalcll2004 


Scalo & Elmegreen 



1999 



2004J). These studies, mainly of our Galaxy and its satellites, find a power law power spec- 
trum P(U) = A U a with index a in the range —2.5 to —3.0. This scale invariant power 
spectrum is interpreted as arisin g from turbulence in the ISM. 

Recently iBegum et al.l ( 120061 ) have presented a visibility based formalism for determining 



the power spectrum of galaxies with ext remely weak HI emission. This 



a sam ple of dwarf and spiral galaxies ( iDutta et al. 



2008 



Dutta et al- 



ias be e n applied to 



2009a : 



Dutta et al. 



2009b|). These studies indicate a dichotomy in a with values ~ —1.5 in some galaxies and 
~ —2.5 in others. The two values ~ —1.5 and ~ —2.5 are interpreted as arising respectively 
from two (2D) and three (3D) dimensional turbulence. Some support for this interpretation 
has been provided by analysis of the fluctuation power spectrum for NGC 1058 where a 
transition from 2D to 3D turbulence is observed at an angular scale corresponding to the 
scale height of the galaxy's disk. 

NGC 4254 the galaxy whose power spectrum we estimate in this paper, is a lopsided 
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Figure 1. The 6.5" X 8.0" resolution integrated HI column density map of NGC 4254 (contours) is overlayed with the optical 
image of the galaxy (density). The contours levels are 3., 5., 10., 20., 25., 35., 40. and 45. xlO 20 atoms cm -2 . 



spira l galaxy (morphological type SA(s)c), with an inclination of ~ 42° ( IPhookun et al. 



19931 ). The galaxy is located at a distance of 1 Mpc from the core of Virgo cluster and is 



believ ed to be falling into the cluster with a relative velocity of 13 00 km s 1 ( Vollmer et al. 



20051 ). The distance to this galaxy is estimated to be 16.7 Mpc (IMei et al. 
distance 1" corresponds to 81 pc. 



20071 ): at this 



2 DATA AND ANALYSIS 

We have used archival HI data of NGC 4254 from the Very Large Array (VLA). The obser- 
vations had been carr ied out on 4 th — 5 th March 1992 using the C configurations of the VLA 



( IPhookun et al.lll993[ ). The data was downloaded from the VLA archive and reduced in the 
usual way using standard tasks in classic aipJ^- The HI emission from NGC 4254 spans over 
30 channels i.e, (18 to 47) from 1401.7 MHz (2253 km s" 1 ) to 1416.4 MHz (2562 km s" 1 ) of 
the 63 channel spectral cube. A frequency width of 48.8 k Hz for each channel in the data 
cube corresponds to the velocity resolution of 10.32 km s _1 . The continuum from the galaxy 
was subtracted from the data in the uv plane using the aips task uvsub. The resulting con- 
tinuum subtracted data was used for the subsequent analysis. Figure [1] shows a integrated 
HI column density (Moment 0) map of NGC 4254 made from this data. The angular ex- 

1 NRAO Astrophysical Image Processing System, a commonly used software for radio data processing. 



4 Prasun Dutta, Ayesha Begum, Somnath Bharadwaj and Jayaram N. Chengalur 

R (kpc) 




tent of the HI distribution in Figured] is measured to be 6.5' x 8.0' at a column density o f 
10 19 atoms cm~ 2 which is co mparable to i t's opt ical diameter (jde Vaucouleurs et al.lll99lf ). 



Begum et al. 



(120061 ) and 



Dutta et al. 



( )2009af ) contains a detailed discussion of the visi- 
bility based HI power spectrum estimator Phi(^) = ( V V (U)V*(U + AU) ), hence we present 
only a brief discussion here. Here U refers to a baseline, i.e. the antenna separation projected 
in the plane perpendicular to the direction of observation, measured in units of the observing 
wavelength A. It is common practice to express the dimensionless quantity U in units of kilo 
wavelength (kA). Every visibility V V (U) is correlated with all other visibilities V*(U + AU) 
within a disk \AU\ < (7t0q) -1 , where 8q is the angular extent of the galaxy. The correlations 
are averaged over different U directions assuming that the signal is statistically isotropic in 
the plane of the galaxy's image. To increase the signal to noise ratio we further average the 
correlation in bins of U and over all frequency channels with HI emission. The expectation 
value of the estimator Pm(U) is real, and it is the convolution of the HI power spectrum 
Pm(U) with a window function |H / V ([/)| 2 . The window function, which quantifies the ef- 
fect of the large-scale HI distribution and the finite angular extent of the galaxy, is peaked 
around U — and has a width of order ( ^0n)~ l beyond which iW^f/)! 2 ~ 0. Analytic con- 



siderations and simulations ( Dutta et al. 



2009al ) show that beyond a baseline U„ 



3.50Q 1 



the convolution does not effect the shape of the HI power spectrum, and we may directly 
interpret the real part of the estimator Pm(U) as the HI power spectrum Pm(U). 

The estimator Pm(U) also has a small imaginary component arising mainly from noise. 
The I — a error-bars for the estimated power spectrum is a sum, in quadrature, of contri- 
butions from two sources of uncertainty. At small U the uncertainty is dominated by the 
sample variance which comes from the fact that we have a finite and limited number of in- 
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dependent estimates of the true power spectrum. At large U, it is dominated by the system 
noise in each visibility. 



3 RESULTS AND DISCUSSION 

Figure [2] shows the real and imaginary parts of Pm(U) evaluated using 16 channels from 
27 to 42 which have relatively high HI emission. Here the estimator Phi(^) was separately 
evaluated for each channel and then averaged over all 16 chann els to increase the s ignal 
to noise ratio. As expected from the theoretical considerations (IBegum et al.l 120061 ). the 
imaginary part is well suppressed compared to the real part. We also estimate Pm({7) using 
the line free channels to test for any contribution from the residual continuum. This is 
found to be much smaller than the signal (Figure [2]) indicating that the continuum has been 
adequately subtracted out. 

The power law Pai(U) = AU a , with a = — 1.7±0.2 is found to give a good fit to the HI 
power spectrum for the U range 2.0 kA to 10.0 kA (Table CD) . The best-fit power- law and the 



1 — a e rror-bar on a were determined using \ 2 minimization ( IBegum et al. 



2006 



Dutta et al. 



2009bl ). We use D — 16.7 Mpc (distance to the galaxy) to convert a baseline U to a length- 



scale D/U in the plane of the galaxy's im age. The U range 2. kA to 10.0 kA corresponds 



to the range of length-scales 8.4 to 1.7 kpc. 



Dutta et al 



(j2009al ) have broadly classified the 



observed turbulence in their galaxy sample as 2D turbulence at the large scales in the plane 
of the galaxies disk and 3D turbulence at scales smaller than the scale height of the disk. 
They also show that the slope of the power spectrum changes from ~ —1.5 to ~ —2.5 as 



is definitely larger than the t ypical HI scale heights within the Mi 



198 



4 



Wouterloot et al. 



ky-Way ( 


Lockman et al. 


Naravan & Joe 


2002 


). It 



19901 ) and in external spiral galaxies (e.g. 
is thus quite reasonable to conclude that the slope a = —1.7 ±0.2 is of 2D turbulence in the 
plane of the galaxy's disk. The fact that we do not o bserve the transitio n to 3D turbulence, 
which is expected to occur at a baseline U = D/nZh (IDutta et al.ll2009al ). allows us to place 
an upper limit on the galaxy's scale height Zh < 2.6 kpc. 

Our present observation is another confirmation (Dutta et al 2008, 2009a 2009b) of the 
fact that spiral galaxies exhibit scale-invariant density fluctuations that extend to length- 
scales of ~ 10 kpc (eg. NGC 628, NGC 1058) which is comparable to the diameter of the HI 
disk. While a large variety of possible energy sources like proto-stellar winds, supernovae, 
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Figure 3. Integrated HI column density maps of the galaxy NGC 4254 using data cubes A, B and C. Note the diagonal 
movement of the centroid of emission from North east (A) to South west (C). The contours levels are 5., 8. and 12. xlO 20 
atoms cm -2 . 



R (kpc) 



iijj A 




NGC 4254 


" f 1 i 








B 


W-JL, 1 ■} 


C 













U(kA) 

Figure 4. HI power spectrum for data cubes A, B and C plotted with arbitrary offsets to prevent them from overlapping. 



shocks, etc. have been proposed to drive turbulence (jEhnegreen fc Scald 12004] ). it is still to 
be seen whether these are effective on length-scales as large as 10 kpc. 

Galaxy harassment is expected to have different effects on the inner and outer parts 
of the galaxy. While gas is stripped from the outer parts, the inner part looses angular 
momentum and gradually collapses to the center through repeated galaxy encounters. We 
can selectively study different parts of NGC 4254, whose rotation axis is tilted at 42° to the 
line of sight, by considering different velocity channels. Our analysis till now has used only 
the central 16 channels, we now use the central 24 channels (23 — 46) for the subsequent 
analysis. We construct 3 different data cubes namely A, B and C containing channels 23 — 30, 
31—38 and 39 — 46 respectively. We can now separately probe the North east, central and 
South west parts of the galaxy (Figure [3]) using these three data cubes. For each data cube, 
we evaluate the HI power spectrum using individual channels and then average over the 
channels in the respective cubes. We estimated U m seperately for A, B and C and the best 
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Figure 5. Integrated HI column density maps of the galaxy NGC 628 using data cubes D, E and F. Note the diagonal movement 
of the centroid of emission from North east (D) to South west (F). The contours levels are 6., 12. and 18. xlO 20 atoms cm . 
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Figure 6. Power spectra of the HI emission for the galaxy NGC 628 is shown for D (channels 108-119), E (channels 120-131) 
and F (channels 132-143) wi t h arb itrary offsets to prevent them from overlapping. Details of the data used for this estimation 
can be found in Dutta et al. (2008) 



fit power law is obtained for U > U m only. We find that for each data cubes the HI power 
spectrum is well fit by a power law (Figure H]), the details being shown in Table [TJ We find 
that the slope a is —2.0 ± 0.3 for A and C which probe the outer parts of the disk while it 
is —1.5 ± 0.2 in B which probes the central region. To verify that this change in slope is due 
to harrass ment, we also consider the HI power spectrum of the spiral galaies NGC 628 and 



NGC 1058 flDutta et al 



2008 



2009bl ) which are not undergoing galaxy harassment. Figure [5] 



shows the regions of the spiral galaxy NGC 628 corresponding to each of the channel range 
108-119 (D), 120-131 (E) and 132-143 (F). We find that the power spectra of these three 
data cubes (Figure [6]) all have the same slope ~ —1.6, which is also similar to the slope of 
the central part of NGC 4254. The results are similar for NGC 1058 and hence we do not 
explicitly show these here. Based on this we conclude that the difference in slope between 
the inner and outer parts of NGC 4254 is a consequence of galaxy harassment. Not only 



8 Prasun Dutta, Ayesha Begum, Somnath Bharadwaj and Jayaram N. Chengalur 



Data 


Channels 


Um 


U range 


a 










(k A) 


(k A) 






16 central 


27 


-42 


1.8 


2.5 - 10.0 


-1.7 ± 


0.2 


channels 














A 


23 


- 30 


1.3 


2.0 - 10.0 


-2.0 ± 


0.3 


13 


31 


- 38 


2.2 


2.5 - 10.0 


-1.5 ± 


0.2 


C 


39 


- 46 


1.3 


2.0 - 10.0 


-2.0 ± 


0.3 



Table 1. Results of the power spectrum analysis of NGC 4254 for different data cubes spanning over different velocity ranges. 

does galaxy harassment affect the global morphology of the galaxy, it also affects the fine 
scale structure in the ISM as reflected by HI power spectrum. 

We note that NGC 4254 have an inclination of 42°, on the other hand NGC 628 and 
NGC 1058 are more face-on galaxies (inclination ~ 10°). In our analysis we selectively study 
different parts of a galaxy by considering different velocity channel ranges. However, since 
face-on galaxies do not have much range in radial velocity from the rotation curve, the spatial 
extent of HI in NGC 628 and NGC 1058 will not be very different for the three velocity 
channel ranges, which is unlike the case for NGC 4254. Hence it is likely that we failed 
to find a change of slope in NGC 628 and NGC 1058 because of this effect of inclination. 
However, as seen in Figure [5]for NGC 628, all three data cubes D, E and F have a significant 
contribution of HI from the outer region of NGC 628, leading us to believe that the power 
spectrum has a slope of ~ 1.6 in the outer parts of NGC 628. This is similar to the value of 
slope seen in the central parts of NGC 4254 and significantly different from the slope in the 
outer parts of NGC 4254. Hence it indicates an impact of harassment in NGC 4254. Further 
analysis of spiral galaxies with large inclination angles would possibly be able to resolve this 
issue. 

We currently do not have an understanding of how galaxy harassment caused a steepening 
of the HI power spectrum in the outer parts of the galaxy. Theoretical modeling and the 
analysis of other Virgo cluster spiral galaxies are needed for further progress in this direction. 
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